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We study the phase structure of SU(3) lattice gauge theory with Nf — 12 staggered fermions 
in the fundamental representation, for both zero and finite temperature at strong gauge couplings. 
For small fermion masses we find two transitions at finite temperature that converge to two well- 
separated bulk phase transitions. The phase between the two transitions appears to be a novel phase. 
We identify order parameters showing that the single-site shift symmetry of staggered fermions is 
spontaneously broken in this phase. We investigate the eigenvalue spectrum of the Dirac operator, 
the static potential and the meson spectrum, which collectively establish that this novel phase is 
confining but chirally symmetric. The phase is bordered by first-order phase transitions, and since 
we find the same phase structure with Nf = 8 fermions, we argue that the novel phase is most likely 
a strong-coupling lattice artifact, the existence of which does not imply IR conformality. 
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I. INTRODUCTION 

Strongly coupled gauge-fermion systems, beyond their 
intrinsic theoretical interest, play an essential role in 
many theories of physics beyond the standard model [1]. 
Lattice gauge theory is at present the most reliable 
method to study strongly-interacting models in a system- 
atic, controlled way. Most lattice studies have focused on 
determining whether given models exhibit confinement 
and chiral symmetry breaking, or if they develop an in- 
frared fixed point (IRFP) resulting in IR conformality [2]. 
SU(3) gauge theory with Nf = 12 fundamental flavors 
has emerged in lattice studies as one of the most contro- 
versial models [3-18]. 

A recent large-scale study of the Nf = 12 system con- 
cluded that the data favored a confining, chirally broken 
scenario [8], though other groups interpreted these data 
as consistent with IR conformality [10, 11]. One of us re- 
cently investigated the same system with an entirely dif- 
ferent approach, the Monte Carlo renormalization group 
(MCRG) two-lattice matching method [6, 9, 15], finding 
an IRFP consistent with IR-conformal dynamics. An ob- 
vious difference between Ref. [8] and Ref. [9] is that the 
MCRG analysis found the IRFP at a considerably weaker 
bare coupling than that considered in Ref. [8]. 

In this work we present a study of the Nf = 12 
system at stronger couplings, reporting results for the 
phase diagram in the gauge coupling-fermion mass pa- 
rameter space, at both zero and finite temperature. We 
find two transitions at finite temperature that converge 
to two well-separated bulk phase transitions, consistent 
with what Refs. [7, 19] observed using different staggered 
lattice actions. The general consistency of results ob- 
tained with very different actions indicates that we are 
observing a robust feature of lattice gauge theories with 
many staggered fermions. Shortly after we completed 
this work, Ref. [14] interpreted the second transition as 
a partial restoration of axial \](1)a symmetry, which is 
not consistent with our data discussed below. 



Between the two bulk transitions we identify a novel 
phase, and will devote most of this paper to understand- 
ing it. We will show that this novel phase breaks the 
single-site shift symmetry of the staggered action, a prop- 
erty that to our knowledge has never been observed be- 
fore. It is important to establish whether or not this 
phase is a lattice artifact, because this will affect the 
conclusions we can draw from our observation of finite- 
temperature transitions converging to bulk transitions. 
As we will discuss in Section III, such behavior is char- 
acteristic of IR-conformal systems, and in principle pro- 
vides a signal that distinguishes confining and conformal 
systems [4]. However, if the phase is a lattice artifact 
such as the "Aoki-like phase for staggered fermions" dis- 
cussed by Ref. [20], then it would be bounded by bulk 
transitions regardless of whether the 12-flavor system 
were confining or conformal in the continuum. More- 
over, some finite-temperature transitions could converge 
to these bulk transitions for both confining and conformal 
systems. As a consequence, our results would not nec- 
essarily imply (although they would be consistent with) 
IR-conformal continuum dynamics. 

The outline and main results of our paper are as fol- 
lows. After summarizing our lattice action in Section II, 
we present our results for the phase structure at light 
fermion mass in Section III. Using two different order 
parameters, in Section IV we will clearly show that the 
single-site shift symmetry ( "S*^" ) of the staggered action 
is spontaneously broken in the novel (".5^") phase. In 
Section V we study the low-lying eigenvalues of the Dirac 
operator in both the ^ phase as well as the more fa- 
miliar weak-coupling phase. The volume scaling of the 
low- lying eigenmodes in the ^ phase indicates the pres- 
ence of a "soft edge" Aq > in the density distribution, 
p(A) (X (A — Aq)". a soft edge implies that not only 
the chiral condensate (V^V^), but also higher-order con- 
densates that could break chiral symmetry, are vanishing 
in the chiral limit [21, 22]. In the weak-coupling phase, 
we do not observe a soft edge in the eigenvalue density 
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distribution. We also obtain a preliminary prediction for 
the mass anomalous dimension 7^ = 0.61(5), where the 
error is purely statistical. 

We present more evidence that the phase is chi- 
rally symmetric in Section VI, through the light meson 
spectrum. In the phase we observe parity doubling be- 
tween scalar and pseudoscalar states as well as between 
vector and axial-vector states. Our meson spectrum re- 
sults in the phase are also independent of the volume, 
in contrast to those in the weak-coupling phase. We also 
investigate the static potential, which shows that the ^ 
phase is confining. The static potential predicts a non- 
vanishing string tension and a small Sommer parameter 
ro ^ 2.7. Minimal finite volume effects both in the meson 
spectrum and static potential are consistent with this ro, 
and show that the lattice correlation length is small in 
^ phase. 

To explore whether the ^ phase could exist without 
implying IR conformality in the continuum, we have per- 
formed preliminary investigations of the phase structure 
with N f = 8 fundamental flavors, a system generally 
believed to show spontaneous chiral symmetry breaking. 
While we have not yet completed our exploration of the 
N j = 8 phase diagram, our initial results in Section VII 
show an phase with the same properties as we ob- 
serve with 12 flavors. We do not observe this phase with 
4 flavors, and our data with 16 flavors are currently too 
preliminary to draw a definite conclusion. Our investi- 
gations with 8 and 16 flavors are ongoing, and we will 
report our full results in future works [23]. 

We conclude in Section VIII. Although we find the 
phase to be confining yet chirally symmetric, it is 
bordered by first-order phase transitions that prevent us 
from reaching the continuum (infinite cut-off) limit. It 
is possible that the B^ phase is a lattice artifact of stag- 
gered fermions, especially since it exists with 8 flavors 
as well. This prevents us from interpreting our observa- 
tion of finite-temperature transitions converging to bulk 
transitions as evidence that the 12-fiavor system is IR- 
conformal. While Ref. [20] has discussed a possible Aoki- 
like lattice artifact phase for staggered fermions, it is not 
clear if the breaking of the single-site shift symmetry we 
observe in the B^ phase can be described by the stag- 
gered chiral Lagrangian. This question deserves further 
study, but it is beyond the scope of the present paper. 



II. THE LATTICE ACTION 

Lattice calculations are affected by discretization er- 
rors, and much effort has been devoted to improving lat- 
tice actions to reduce these effects. Strongly-coupled sys- 
tems are particularly sensitive to these lattice artifacts, 
which can contaminate or destroy the scaling of the de- 
sired continuum limit, even to the point of generating 
spurious ultraviolet fixed points (UVFPs). Care must 
be taken that lattice simulations are in the basin of at- 
traction of the perturbative fixed point, or its associated 



IRFP if it exists. 

In Ref. [9] we advocated the use of a gauge action with 
a negative adjoint plaquette term, to avoid a well-known 
spurious UVFP caused by lattice artifacts. In the present 
work we follow this suggestion and use a gauge action 
that includes both fundamental and adjoint plaquette 
terms, with coefficients related by = —0.25/3^?. With 
this constant ratio, the perturbative relation to the bare 
gauge coupling is 

6/^2 = + 2/3a//3^) = /3i./2. (1) 

Staggered fermions are also affected by taste break- 
ing, i.e., the four fermion tastes described by each (un- 
rooted) staggered field are degenerate only in the con- 
tinuum limit. Smearing the gauge connections reduces 
this problem, and following Refs. [6, 9, 24], we use 
nHYP-smeared staggered fermions. The nHYP smear- 
ing significantly improves the taste symmetry of stag- 
gered fermions [25, 26], but the U(3) projection in the 
nHYP construction can break down at strong coupling, 
due to the generation of near-zero eigenvalues in the sta- 
ple sum. We address this difficulty by adjusting the three 
HYP smearing parameters to 

Oi\ = 0.5, OL2 = 0.5, as = 0.4 

from the original (0.75, 0.6, 0.3) values. This choice elim- 
inates the numerical problems while it only slightly in- 
creases taste splitting. 

In our calculations we use the hybrid Monte Carlo 
(HMC) algorithm. Our code is based in part on the 
MILC Collaboration's public lattice gauge theory soft- 
ware [27]. We have modified this software to implement 
nHYP smearing, to add the adjoint plaquette term to 
the gauge action, and to exploit both even and odd sub- 
lattices to simulate eight flavors. During the course of 
our work, we also implemented a second-order Omelyan 
integrator [28] accelerated by an additional heavy pseud- 
ofermion field [29] and multiple time scales [30]. Our 
HMC trajectory length is typically one molecular dynam- 
ics time unit (MDTU), but in some cases can be as small 
as 0.5 MDTU or as large as 2.0 MDTU. For most of the 
ensembles discussed in this work, we accumulated 1000- 
2000 MDTU, and measured the eigenvalues and meson 
spectrum on every tenth trajectory. Around the phase 
transitions we accumulated more than 10,000 MDTU for 
some ensembles. 



III. THE PHASE STRUCTURE 

In the m = chiral limit, confining and chirally broken 
systems with Nf > 3 flavors of fundamental fermions 
are expected to exhibit a first-order finite-temperature 
phase transition at which they become chirally symmetric 
and deconfined. A finite-temperature lattice system with 
fixed Nf <^ L will undergo a phase transition at a critical 

(c) 

coupling . In the weak-coupling scaling region the 
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FIG. 1. The chiral condensate (ipip) (on a log scale) and 
the blocked Polyakov loop {Ti Lb) indicate two well-separated 
transitions at m = 0.005. 



renormalization group equation predicts the dependence 
of 13 on Nf. In order for the theory to be confining 

and chirally broken at zero temperature, (3^"^ ^ oo as 
Nt oo. 

Zero-temperature systems with Nt > L deconfine and 
become chirally symmetric when L is so small that the 
physics is volume-squeezed. This is a finite-volume ef- 
fect and not a real phase transition, though it could be 
accompanied by a discontinuity. In few-fiavor QCD-like 
systems no discontinuity is observed at zero temperature. 

Much less is known about the finite-temperature be- 
havior of IR-conformal systems. At strong enough cou- 
pling, lattice artifacts create a confining, chirally broken 
phase on the lattice [31]. This strong-coupling phase 
must be separated from the weak-coupling conformal 
phase by a bulk (non-thermal) phase transition in the 
chiral limit. The bulk transition has to be a real infinite- 
volume transition, with the chiral condensate (V^V^) serv- 
ing as an order parameter in the chiral limit. While 
remnants of the finite-temperature phase transition can 
coexist with the bulk transition, the finite-temperature 
transitions must occur at stronger couplings than the 
bulk transition, and converge to the bulk transition as 
Nt ^ oo. This, in principle, gives a signal that distin- 
guishes confining and conformal systems. 

Refs. [4, 7, 14] investigated the 12-fiavor SU(3) system 
and found indication for a bulk transition. Refs. [7, 14, 
19] also discussed a second discontinuity in ('0'^/^). Our 
investigations confirm the existence of both bulk phase 
transitions, as illustrated in Fig. 1. In the chiral conden- 
sate we observe a clear discontinuity around Pf ^ 2.0 for 
m = 0.005 on zero-temperature 8^, 12^ and 16^ volumes. 
{iptlj) has another, much smaller, discontinuity around 
$F ^ 2.65, where the Polyakov loop, an observable re- 
lated to confinement, shows a much stronger signal. Be- 
cause the usual Polyakov loop becomes small and noisy 
as Nt increases, we consider an improved observable by 
measuring the Polyakov loop on renormalization group 
blocked lattices. This blocked Polyakov loop (TrL^) has 



the same Zs symmetry as the standard one, and can also 
be thought of as an extended observable on the original, 
unblocked lattices. 

The chiral condensate (V^V^) has very little volume de- 
pendence across the phase transitions, consistent with 
bulk transitions. The apparent volume dependence of 
the blocked Polyakov loop is due to the different number 
of blocking steps performed: the 16^ lattices are blocked 
three times with scale factor s = 2, while the 12^ lattices 
can be blocked only twice. 

Fig. 2 shows how the locations of the two phase tran- 
sitions depend on the volume, temperature and fermion 
mass. In both cases, the transitions at finite tempera- 
ture converge to zero-temperature bulk transitions where 
different observables show the same discontinuity on all 
volumes, up to small finite volume effects. Just like 
Refs. [4, 7, 14], we observe the stronger-coupling tran- 
sitions to converge on smaller volumes than those that 
are needed for the weaker-coupling transitions to con- 
verge. We encountered long metast ability between runs 
from hot and cold initial states at both transitions, espe- 
cially on larger volumes. These transitions are strongly 
first-order, and molecular dynamics evolution is not very 
effective fiipping the system between phases. Mixed ini- 
tial configurations helped to resolve the transition around 
Pf ^ 2.65 more accurately, but they were less reliable at 
the transition around /3f ~ 2. 



IV. SINGLE-SITE SHIFT SYMMETRY 
BREAKING 



We identified two phase transitions in Fig. 1 from the 
discontinuities in the chiral condensate and (blocked) 
Polyakov loop. This was possible as both phase tran- 
sitions are first-order, and almost all observables show a 
discontinuity. However, neither (V^V^) nor the Polyakov 
loop is a bona fide order parameter of the intermedi- 
ate phase located between the two transitions. While 
the Polyakov loop is a good indicator of confinement, 
it is only an order parameter in the pure gauge theory, 
and does not distinguish between the intermediate and 
strong-coupling phases. The chiral condensate is an or- 
der parameter in the chiral limit only, and in that limit 
it likely vanishes in both the intermediate and weak- 
coupling phases. 

There is no a priori guarantee that the intermediate 
phase is separated from the strong- and weak-coupling 
phases by true phase transitions. However, while investi- 
gating the phase diagram, we discovered that the single- 
site shift symmetry C'^^") of the staggered fermion ac- 
tion is spontaneously broken in the intermediate ( " ) 
phase. This ensures the existence of an order parameter 
characterizing the J8^ phase, and full separation of the 
phases. 

The single-site shift symmetry of the staggered action 
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FIG. 2. Transitions in the Pf-tti plane at several temperatures and volumes. The left panel shows the stronger-coupling 
transitions signaled by (ipip), the right panel shows the weaker-coupling transitions signaled by (TrL^), and the middle panel 
shows how the two bulk transitions merge as m increases. Small vertical offsets distinguish the different volumes, and lines 
connect the points to guide the eye. The transitions are nearly identical on zero-temperature volumes, and the finite-temperature 
transitions appear to converge to these bulk transitions. 
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FIG. 3. The plaquette difference APt (Eqn. 3, left) and the link difference ALt (Eqn. 4, right) measured on 16x32 volumes 
in both the phase [Pf — 2.6, m — 0.005) and the weak-coupling phase [I^f — 2.7, m — 0.005), as functions of the molecular 
dynamics time. 



takes the form [32] 



where 



(2) 



C^(n) = (-1)^-.' 



This symmetry ensures that the chiral condensate (V^V^) 
measured on even lattice sites is identical to that mea- 
sured on odd sites, and the underlying gauge configu- 
rations exhibit the usual discrete translational symme- 
try. To our knowledge the breaking of this symmetry has 
never been observed before. 

Order parameters that are sensitive to this symme- 
try include the expectation value of the difference be- 
tween neighboring plaquettes and that between neigh- 
boring links, 



AP^ = (Re Tr - Re Tr □n+M)n, 
AL^ =(a^(n)x(n)/7^(n)x(n + /i) 



(3) 
(4) 



+ m)x(^ + l^)Uf,{n + /i)x(n + 2/i))n^e 



where Re Tr is the real trace of the plaquette origi- 
nating at site n, 

is the usual staggered phase factor, and the expectation 
value (• • even taken only over sites whose /i compo- 
nent is even. In the intermediate phase these operators 
develop non-zero expectation values in one or more di- 
rections /i. Occasionally the direction of the symmetry 
breaking changes, rotating in space. 

Fig. 3 shows both order parameters in the intermediate 
{fSp = 2.6, m — 0.005) and the weak-coupling = 2.7, 
m = 0.005) phases, as functions of the molecular dynam- 
ics time on 16^x32 volumes. At = 2.6 is broken in 
the temporal direction and both order parameters remain 
small for fi ^ t. At fSp = 2.7 the order parameters do 
not develop a non-zero expectation value in any direction. 
The order parameters also vanish in the strong-coupling 
phase, though we do not include that data in Fig. 3. The 
single-site shift symmetry is broken only in the interme- 
diate phase. The order parameters, when non- vanishing, 
have only small dependence on the volume. 
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It is important to note that the single-site shift sym- 
metry is an exact symmetry of the action even at finite 
fermion mass. It is broken only spontaneously. Both 
AP^ of Eqn. 3 and AL^ of Eqn. 4 are nonzero when 
the symmetry is broken and vanish when it is preserved. 
The JS^ phase must be separated by true phase transi- 
tions from the strong- and weak-coupling phases where 
both order parameters vanish. 



V. EIGENVALUE SPECTRUM 



where the parameter Aq > allows the possibility of a 
soft edge [21, 35, 36]. In a chirally broken system p(0) / 
0, implying Aq = (the "hard edge") and a = 0. In a 
conformal system Aq = 0, and a is related to the mass 
anomalous dimension as we now derive. 

Although the density p{X) is only well-defined in the 
infinite- volume limit, the functional form of Eqn. 5 can be 
used to analyze the spacing between discrete eigenvalues 
in a finite volume. Denoting the finite- volume eigenval- 
ues as Ai for i = 1, 2, . . . , we can write the cumulative 
eigenvalue density as 



The results we discussed in Section III suggest that the 
transition at stronger coupling is related to chiral sym- 
metry breaking, while the transition at weaker coupling 
is related to confinement. We will consider confinement 
in the next section, while in this section we investigate 
the chiral properties of the phases. Finite fermion mass 
explicitly breaks chiral symmetry, and extrapolating the 
chiral condensate to the m = chiral limit can be a diffi- 
cult task. Here we use the eigenvalue distribution of the 
Dirac operator to study chiral symmetry. 

The spectrum of the Dirac operator of chirally bro- 
ken systems contains a wealth of information. When the 
eigenvalue distribution is compared to random matrix 
theory (RMT), it predicts the chiral condensate and gives 
information about the lattice artifacts of the simulations. 
There is no comparable prediction for conformal systems, 
nevertheless the volume scaling and the level spacing of 
consecutive eigenvalues are related to the mass scaling 
exponent of the fixed point that governs the infrared dy- 
namics [33, 34]. 

The 12-fiavor staggered action is local and describes 
a well-defined statistical system. In the following anal- 
ysis we will investigate general questions of scaling and 
the appropriate scaling dimension for this system. We 
will not compare our data to RMT predictions, and our 
analysis will not be sensitive to taste symmetry, or its 
breaking. 

In this pilot study we calculate the 12 lowest-lying 
eigenvalues of the staggered Dirac operator on volumes 
12"^, 12^x24, 16"^ and 16^x32 in both the ^ and weak- 
coupling phases. In principle one should separate the dif- 
ferent topological sectors before averaging the eigenval- 
ues, but all of the configurations we analyzed appear to 
be in the zero-topology sector. Fig. 4 illustrates the vol- 
ume dependence and the level spacing of the lowest-lying 
eigenvalues. Our gauge configurations are too coarse for 
the eigenvalues to show the four-fold degeneracy expected 
in the continuum limit of staggered fermions. Additional 
HYP smearing steps can remove enough of the ultravi- 
olet fiuctuations to reveal the expected degeneracy, but 
that is not what the dynamical fermions see in the sim- 
ulations, and we do not pursue this direction. 

In the infinite-volume limit the basic quantity is the 
eigenvalue density p(A). In the chiral limit the density of 
low-lying eigenvalues is expected to scale as 



p{X)dX 



X 



lim 



V 



(6) 



where A^ = A and A^ = A. Using Eqn. 5 this leads to 

1 

n — xq ^ 



An - Ao (X 



V 



[l + (9(F-i)], (7) 



where we combined m/V and {Xm — Aq)"^"^ cx 0{V~^) 
in the parameter xq/V. Since the eigenvalue A has di- 
mension of mass, Eqn. 7 implies the relation 



D 



1 



(8) 



between a and the mass anomalous dimension in a D- 
dimensional space. For free field theory = 1 and a = 
D — 1^ while for a chirally broken system a = 0. Eqn. 7 
has four free parameters. The proportionality constant 
and xq can depend on the lattice geometry, while Aq and 
a are universal. 

The left panel of Fig. 4 shows our results for the low- 
lying eigenvalues in the phase at/3F = 2.6,m = 0.005. 
In this plot we include a dashed line showing the soft edge 
Ao = 0.0175(5) predicted by our global fit to Eqn. 7 using 
all four volumes. The left panel of Fig. 5 shows the result 
of this global fit, which does not depend on the aspect 
ratios of the lattices. The dependence on xq is weak, and 
we fix xo 0. The slope ym 4/(a + 1) = 2.50(10) 
gives a = 0.60(6), consistent with the RMT prediction 
a = 1/2 at a soft edge [35, 36]. 

A non- vanishing soft edge is rather unusual. In finite- 
temperature systems with Nf fixed, L ^ 00, Ref. [21] ob- 
served Ao > in the chirally broken phase, but in infinite 
volume neither chirally broken nor conformal systems are 
expected to have a soft edge. Through the Banks-Casher 
relation [37] 



Jo 



p{X)dX 
A^ 



(9) 



(5) 



a soft edge implies that the chiral condensate (V^V^) van- 
ishes in the chiral limit m = 0. With a soft edge, p(A) = 
for < A < Ao as well as for all A larger than the spectral 
range of the Dirac operator, so that the integral in Eqn. 9 
remains finite while m ^ 0. 
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FIG. 4. The low-lying eigenvalues (A^) for m = 0.005 on the four volumes 12^, 12^x24, 16^ and 16^x32, in 
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FIG. 5. Scahng of the low- lying eigenvalues for m = 0.005 on the four volumes 12^, 12^ x 24, 16^ and 16^ x 32. In the ^ 
phase = 2.6, left) the soft edge Ao = 0.0175(5), and the slope ym = 2.50(10). In the weak-coupling phase (/3f = 2.7, right), 
Ao = but xo ~ 3, with slope ym — 1.61(5). 



In addition, a soft edge excludes the scenario in which 
(ipip") = but chiral symmetry is broken in the phase 
by a nonzero four-fermion condensate. As discussed by 
Refs. [38-40], this could result from the chiral symmetry 
breaking pattern 

SU{Nf)v X SU{Nf)A ^ SU{Nf)v x Zn^ (10) 

where the custodial Z^f symmetry forces {^I^iIj) = 0. The 
four-fermion condensate considered in Ref. [40] is related 
to the difference of scalar and pseudoscalar susceptibili- 
ties uj = xp — Xs where 



Xp 



Xs 



^ J d^xd^y{xrH-f5x{x)xr^n5X{y)) (H) 
= ^J d'xd'y{xr^x{x)xr'x{y)) (12) 



and is a flavor generator. The U(l)A-noninvariant uj 
parameter can be expressed in terms of the eigenvalue 
density as [22] 



u 4m^ 



f 

Jo 



p{X)dX 
(A2 + m2)2 



(13) 



Just as for Eqn. 9, uj vanishes in the chiral limit if 
the eigenvalue density has a soft edge, so the symme- 
try breaking scenario of Eqn. 10 is not consistent with 
our data in the phase. 

The soft edge is a dimensional parameter, but it is 
not clear what infinite- volume physical quantity it corre- 
sponds to. Better understanding of the symmetry break- 
ing mechanism could shed light on this question. 

The eigenvalue spectrum in the weak-coupling phase 
is more conventional. The right panel of Fig. 4 shows 
the low- lying eigenvalues in this phase at = 2.7, 
m = 0.005. The global fit to Eqn. 7 predicts Aq = but a 
non- vanishing xq ^ 3. Volumes with different aspect ra- 
tios prefer slightly different Xq values and proportionality 
constants. A fit with a common a parameter to the 12 
eigenvalues on all four volumes predicts y^ = 1.61(5) or 
mass anomalous dimension 7^ = 0.61(5) (right panel of 
Fig. 5), where the error is purely statistical. This value 
is consistent with results reported by Refs. [10, 11, 17]. 
We find similar scaling properties and exponent at mass 
m = 0.01. 

In a chirally broken, confining system the eigenvalues 
should scale with exponent ym = 4, even in the e regime 
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The meson spectrum at = 2.7 is also consistent 
with a small- volume deconfined system. The right panel 
of Fig. 7 shows the Goldstone 75 pseudoscalar (tts), the 
pseudoscalar and the scalar components of the 7075 chan- 
nel (ttos and ao) and the 7^75 pseudoscalar (tt^s) versus 
fermion mass m. We observe significant volume depen- 
dence in the scalar and pseudoscalar masses. The Gold- 
stone and the scalar become degenerate at small m, con- 
sistent with parity doubling. The ttqs meson becomes 
heavier than the scalar at m = 0.005, where it is degen- 
erate with the 7r^5 state. (Our data do not allow precise 
results for these states on 12^x24 at m < 0.01.) Overall 
our meson spectrum results at = 2.7 are dominated by 
finite- volume effects, and do not provide clear informa- 
tion about the IR dynamics of the 12-fiavor model. With 
the computational resources available to us, we cannot 
compete with the large- volume spectral study of Ref. [8] . 



where the volume is small compared to the pion Comp- 
ton wavelength. Our observation of uniform scaling on 
all four volumes indicates that the 12-fiavor system is ei- 
ther conformal, or that its intrinsic confinement scale can 
only be observed on larger lattice volumes than we con- 
sider here. In future work we will obtain a more robust 
prediction and quantify systematic effects by performing 
similar calculations at different gauge couplings and mass 
values, and on larger volumes as well [23]. For now, we 
turn to investigating the confinement properties of the 
and weak-coupling phases. 



VI. STATIC POTENTIAL AND MESON 
SPECTRUM 

In this section we explore the static potential and me- 
son spectrum in the JS^ phase, contrasting these results 
with the same observables in the weak-coupling phase. 
Although (TrL^) shows a clear signal in Fig. 1, the 
Polyakov loop is not an order parameter in the presence 
of dynamical fermions. The static potential is a more 
reliable indicator of confinement. In Fig. 6 we contrast 
the HYP-smeared static potential [41] measured on each 
side of the transition, at I3f = 2.6 and 2.7 on 12^x24 and 
16^ X 32 volumes with m = 0.005. 

The potential at Pf = 2.6 is consistent with confine- 
ment, with string tension a = 0.20(1) and Sommer pa- 
rameter ro = 2.67(4) in lattice units. We obtain similar 
values at other masses and couplings within the >S^ phase. 
The potential is almost identical on 12^x24 and 16^x32 
volumes, and the small ro suggests that there will be no 
qualitative change on larger volumes that we are cur- 
rently investigating. These results indicate confinement 
with a fairly short gauge correlation length. On the other 
hand, the potential at = 2.7 is coulombic and cannot 
be fitted consistently with a linear term. The lack of vol- 
ume dependence implies either vanishing string tension 
and conformality or an intrinsic confinement scale that 
can only be observed on larger lattice volumes. 



Our goal in investigating the static potential and me- 
son spectrum in the weak-coupling phase is to contrast 
these results with measurements in the JS^ phase, where 
we observe several interesting differences. Our results for 
the pseudoscalar and scalar spectrum at /3f = 2.6 are 
summarized in the left panel of Fig. 7. In the phase 
we find that the pion has a parity partner ("as") in the 75 
channel, a state that is forbidden in QCD-like systems. 
The masses measured on 16^ x 32 and 12^ x 24 volumes 
are indistinguishable in the JS^ phase: the finite volume 
corrections are negligible, consistent with the small cor- 
relation length indicated by the static potential. The 
parity partner states both in the 75 and 7075 channels 
are degenerate. The 75 states are largely independent 
of the fermion mass m while the 7075 mesons' masses in- 
crease steadily with increasing m. The data indicate that 
all four mesons could be degenerate in the chiral limit. 
However, the ttqs mass again proved difficult to extract, 
and our statistics and volumes do not allow precise re- 
sults for the TTos at m > 0.01. 

In Fig. 8 we show the masses of the vector meson p and 
its parity partner ai measured on 16^x32 volumes in both 
the phase at (3f = 2.6 and the weak-coupling phase 
at (Sf = 2.7. At both of these couplings, the p and ai are 
degenerate for all m < 0.015. In the deconfined weak- 
coupling phase, this parity doubling is a familiar effect. 
In the confining phase, however, such spectral proper- 
ties are unusual. The p-ai parity doubling we observe in 
Fig. 8 is inconsistent with the spectrum associated with 
the chiral symmetry breaking pattern of Eqn. 10 [40]. 
Combined with the vanishing chiral condensate observed 
from the Dirac spectrum in this phase, the degeneracy of 
the parity partners in the meson spectrum implies that 
the intermediate phase is confining but chirally symmet- 
ric. The continuum 't Hooft anomaly matching condition 
does not permit the existence of such a phase, suggest- 
ing that the novel phase we observe does not exist in the 
continuum. 
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FIG. 7. The masses of light scalar and pseudoscalar staggered mesons, from 12^ x 24 and 16^ x32 lattices. Small horizontal 
offsets distinguish results from different volumes. In the left panel we include the Goldstone tts, its "as" parity partner, the 
TTos pseudoscalar and the ao scalar in the intermediate phase at /3f = 2.6. In the right panel we show the tts, ttos, ao and (on 
16^ x32 only) in the weak-coupling phase at pF = 2.7. 
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FIG. 8. The p and ai are degenerate in both the weak- 
coupling phase at 13 f — 2.7 as well as the phase at /3f = 
2.6. At pF = 2.7 this degeneracy is a familiar effect. At /3f = 
2.6, it is consistent with our observation of chiral symmetry 
in the eigenvalue spectrum. 



FIG. 9. The chiral condensate (on a log scale) in the 

Nf — 8 flavor system at m = 0.005 on 12^^ and 16^ lattices. 
The phase between the two first order transitions is an 
phase like that we observe for Nf — 12. 



VII. THE 8-FLAVOR CASE 

Finite-temperature transitions converging to a bulk 
transition could signal that the continuum weak-coupling 
phase is conformal in the infrared. However, because we 
observe two bulk transitions bounding an intermediate 
phase with unusual properties, we must consider the pos- 
sibility that our results are due to lattice artifacts. With 
Wilson fermions the existence of a lattice artifact phase, 
first proposed by Aoki [42], is well known. Ref. [20] ar- 
gues that an Aoki-like phase might exist with staggered 
fermions if more than a single four-taste multiplet is con- 
sidered. We are currently investigating this possibility 
through additional studies with Nf = 4, 8 and 16 fla- 
vors. 

This work is preliminary, but important in interpret- 
ing our Nf = 12 results. Fig. 9 shows the Nf = 8 chiral 
condensate (V^V^) at m = 0.005 on 12^ and 16^ volumes. 



We observe the same phases as with Nf = 12 flavors. On 
both volumes there are two first-order transitions at ap- 
proximately volume-independent gauge couplings. The 
phase in between has the same properties as the phase 
with 12 flavors. It breaks single-site shift symmetry as 
shown by the non-zero expectation values of the two or- 
der parameters AP^ (Eqn. 3) and AL^ (Eqn. 4). The 
Dirac operator eigenvalue spectrum has a soft edge, the 
static potential has a non-vanishing string tension, and 
the meson spectrum shows parity doubling. Yet it is gen- 
erally believed that the Nf = 8 flavor system is below the 
conformal window [3, 5, 33, 43, 44], and our data in the 
weak-coupling phase support this expectation [23]. 

The fact that an JS^ phase exists with 8 flavors implies 
that this phase and its two corresponding bulk transitions 
do not necessarily imply IR conformality in the contin- 
uum. The infrared behavior of the weak-coupling phase 
is independent of the phase and has to be studied by 
other means. 
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We emphasize that our treatment of the Nf = 8 system 
is not complete, and our results with Nj = 4 and 16 
flavors are even more preliminary. We include Fig. 9 to 
help clarify the situation with 12 flavors, but more work 
is needed to map out the full phase diagram. 



VIII. CONCLUSION 

Our investigations of the phase diagram of the 12-flavor 
SU(3) model have identifled a novel phase with unusual 
properties. At small masses this phase lies in between the 
usual conflning, chirally broken lattice strong-coupling 
phase and the weak-coupling phase that is governed by 
the perturbative gaussian flxed point and possibly an in- 
frared flxed point. The two flrst-order phase transitions 
separating these three phases get closer together with in- 
creasing fermion mass. At some mass value the two tran- 
sitions merge and eventually turn into a crossover. The 
intermediate phase forms a packet in between the strong- 
and weak-coupling phases at small fermion masses. 

In this work we studied the intermediate phase and 
contrasted it with the weak-coupling phase using several 
observables. The chiral condensate (V^V^) and blocked 
Polyakov loop gave our flrst glimpse of the phase struc- 
ture, and suggested that the transition at stronger cou- 
pling is related to chiral symmetry breaking, while the 
transition at weaker coupling is related to conflnement. 
Our investigation led us to two operators, AP^ (Eqn. 3) 
and AL^ (Eqn. 4), that serve as order parameters of 
the intermediate phase. Both of these order parame- 
ters are sensitive to the single-site shift symmetry (S^) of 
the staggered fermions, a symmetry that is exact at the 
level of the lattice action. Since these order parameters 
develop non-zero expectation values in the intermediate 
phase, but vanish in both the strong- and weak-coupling 
phases, we conclude that the intermediate phase sponta- 
neously breaks single-site shift symmetry, JS^. Since the 
single-site shift symmetry is exact even at flnite fermion 
mass, the phase must be separated by real phase tran- 
sitions from both the strong- and weak-coupling phases. 

We used the spectrum of the Dirac operator to study 
the chiral properties of the phases. In the JS^ phase we 
found evidence for a soft edge, implying chiral symme- 
try. In the weak-coupling phase the eigenvalue spectrum 
is consistent with both conformal and volume-squeezed 
conflning scenarios. We obtained a preliminary predic- 
tion for the mass anomalous dimension, 7^ = 0.61(5) in 
the weak-coupling phase. 

The static potential showed that the JS^ phase is con- 
flning with a small lattice correlation length, while in 



the weak-coupling phase on our relatively small volumes 
the potential was only coulombic. These results are con- 
sistent with the signal from the (blocked) Polyakov loop. 
The meson spectrum in both phases indicated parity dou- 
bling at light fermion mass. However, in the JS^ phase 
we observed very little volume dependence, yet we found 
that all mesons remained massive in the chiral limit. The 
parity doubling in the weak coupling phase was accom- 
panied by strong volume dependence and could also be 
consistent with both conformal and conflning scenarios. 

We presented preliminary data showing that the JS^ 
phase is present with 8 flavors as well, suggesting that 
this phase is not related to conformal infrared dynamics. 
Our flndings lead us to believe that the phase is a 
lattice artifact of the staggered fermions [20]. Since the 
single-site shift symmetry is closely related to the fermion 
staggering and taste breaking, it is most likely that the 
origm ofthe^ph ase is in the fermionic sector. However, 
we do not yet have a clear and complete understanding 
of the symmetry breaking mechanism that produces the 
phase. 

Further investigations of the phases on larger volumes, 
and with Nf = 4, 8 and 16 fermions, are under way and 
should clear up the still open questions of these surpris- 
ingly complex systems. 
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